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Abstract 

This article reports the optimized experimental requirements to determine neutrino mass hier¬ 
archy using electron antineutrinos (Pe) generated in a nuclear reactor. The features of the neutrino 
mass hierarchy can be extracted from the and |Am 32 | oscillations by applying the Fourier 

sine and cosine transform to the L/E spectrum. To determine the neutrino mass hierarchy above 
90% probability, the requirements on the energy resolution as a function of the baseline are studied 
at sin^ 2013 = 0.1. If the energy resolution of the neutrino detector is less than Ep and the 

determination probability obtained from Bayes’ theorem is above 90%, the detector needs to be 
located around 48-53 km from the reactor(s) to measure the energy spectrum of i>e- These results 
will be helpful for setting up an experiment to determine the neutrino mass hierarchy, which is an 
important problem in neutrino physics. 
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I. INTRODUCTION 


Since the measurement of the large sin^ 26^13 at RENO, Daya Bay, and Double Chooz, 


the precise measurement o 
focus in neutrino physics [l- 


neutrino mass hierarchy, the sign of Am32, has become the 
3 |. It had been believed that the neutrino mass hierarchy can 
be determined through long-baseline experiments, mainly using accelerator neutrino beams. 
Recently, the capability of a reactor neutrino experiment at an intermediate baseline to 
distinguish normal or inverted hierarchy was reported. 

For an intermediate-baseline neutrino experiment, many approaches have been proposed; 


they can be categorized into the analysis methods, which are discussed in Refs 


and the Fourier-transform methods 


■Ml, 


10l |. The analysis methods based on the newly 


adopted Bayesian approach utilizes all the available information from experiments, and 
it is straightforward to incorporate the uncertainties in order to evaluate the sensitivity, 
providing robust and complementary results in the Fourier-transform methods [ h]. Although 
the analysis methods are attractive and interesting, the Fourier-transform methods are 
more intuitive. The prominent merit of the Fourier-transform methods are that the mass 
hierarchy can be extracted without precise knowledge of the reactor antineutrino spectrum, 
the absolute value of the large [Amg^l, and the energy scale of a detector. The Fourier- 
transform methods were introduced to enhance and visualize the structures of mass hierarchy 
in the frequency spectrum, as hrst discussed in Ref. |l 2 |. 


In principle, the mass hierarchy can be determined through precise measurements of 
|Am3;^| and |Am32|. As lAm^^l is very small and is only ~ 3 % of [Amg^^l, we have to 
measure lAmg^l and |Am|2| with a precision much better than 3 %. However, |Am|]^| and 
|Am32| have been measured in many experiments with only S> 3 % precision (l 3 |. 

The intermediate baseline based on reactor neutrino experiments has been explored using 
the precise measurement of distortions of the energy spectrum with negligible matter effect. 
Learned et al. proposed a new method to distinguish normal and inverse hierarchy after 
a Fourier transform of the L/E spectrum of reactor neutrinos 12 j . They pointed out that 
the Fourier power spectrum has a small but not negligible shoulder next to the main peak, 
and its relative position could be used to extract the mass hierarchy while a non-zero 6*13 is 
considered. 

In this paper, we analyze the sensitivity of medium-baseline reactor antineutrino exper- 
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iments to the neutrino mass hierarchy for a baseline range of 30-60 km and overall energy 
resolution, 5E/ in the range of 0 to 0.08/\/^ with the Fourier-transform method. The 
optimal baseline length is estimated based on the expected probability of determination. 


II. DETECTION OF REACTOR ANTINEUTRINO 


In a nuclear reactor, antineutrinos are mainly produced via the /9-decay of the hssion 
products of the four types of radioactive isotopes, 235^^^238 [ 7^239 in the fuel. 

The antineutrino flux having energy Ey in MeV with thermal power P^h in GWth is given as 


dN 


Pi 


th 


-(j){Ey) X 6.24 X 


( 1 ) 


dEy fk^k 

where fk and ek are the relative hssion contribution and the energy released per hssion of 
isotope k, respectively. Further, 1 is the number of antineutrinos produced per hssion 
and is obtained as follows 1% Il6l |: 


4>{Ey) = /235[/e' 


0.870-0. 16 OE 1 .-0.091^2 






0.793-0.080E^-0.1085EH 


( 2 ) 


The antineutrino hux is modulated by neutrino oscillation. The antineutrino survival prob¬ 
ability Pee is expressed as 


Pee = 1 - cos^( 6 'i 3 ) sin^(26'i2) sin^(A 2 i) 

- cos^(6'i 2) sin^(26'i3) sin^(A3i) 

- sin2(6'i2) sin2(26'i3) sin2(A32). (3) 


The oscillation phase Aj,- is dehned as 


Ajj — 


ArrifjL 
iEy ’ 


/ A 2 _ 2 2\ 

[Ami. =mi- mA 


( 4 ) 


with a baseline L. As A 31 and A 32 appear simultaneously in Eq. (3), the ehects of mass 
hierarchy on Pee are hardly recognized. By using the relation between the squared mass 
diherences, 

= 0 ; (5) 
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we rearrange Eq. (3) to eliminate the A 32 term as follows: 


Pee = 1 - cos‘^(6'i3) sin^(26'i2) sin^(A 2 i) 

- sin^( 26 'i 3 )sin 2 (A 3 i) 

- sin^(6'i2) sin^(26'i3) sin^(A2i) cos(2|A3i|) 

± ^^^^^sin^(26»i3)sin(2A2i)sin(2|A3i|). (6) 

The plus (minus) sign in the hfth term on the right-hand side of Eq. ( 6 ) corresponds to the 
normal (inverted) mass hierarchy or NH (IH) in short. 

In ongoing reactor experiments, we assume that protons will be used as targets to de¬ 
tect electron antineutrinos via the inverse-beta-decay (IBD), which produces a neutron and 
positron. The antineutrino distribution observed with a detector having Np free protons can 
be expressed for an exposure time T as follows: 


dEy 


NpT dN 

47rL2 dE^ 


Pee^L, Ey^CTippi^Ep) ^ 


(7) 


where (Tibd is fhe cross section of the IBD process and L is the baseline length. 

We use the distribution of the expected antineutrino events from the above expression. 
For the IBD cross section, we use the following expression from Vogel and Beacom’s workjl^: 


(^iBD = 0.0952(PePe) X 10 "^^cm^. 


( 8 ) 


In order to study the sensitivity of the mass hierarchy, we use Fourier-transform method 
together with Monte-Carlo simulations to compare the simulated IBD energy spectrum with 
the expected spectrum in both the NH and IH cases. 

Taking into account the detector response, the reactor electron antineutrino z/g L/E 
spectrum becomes 


dE°^^ 



dN°^^ 

dE^ 


R{E^, E^ 


obs 

V '! 


SE,), 


(9) 


where E^ is the actual z/g energy, E°^^ is the observed z/g energy with the detector response, 
5Ey is the energy resolution, and P(P, E') describes the detector response function including 
effects such as the energy resolution and energy scale. In this study, we take the normalized 
Gaussian function as the response function: 


R{E°J>%E,,6E,) = 


\/^5Ei 


-exp ■ 


(E. 


obs 

V 




25E„ 


( 10 ) 
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As the neutrino energy is usually measured using scintillators, the energy is typically pro¬ 
portional to the number of photoelectrons, and the error is dominated by the photoelectron 
statistics. Therefore, the neutrino energy resolution is proportional to In general, 

the detector energy resolution is parameterized into two parts: 


6E, 

E, 



( 11 ) 


The hrst term represents the uncertainty from statistical fluctuation, and the second term 
originates from the systematic uncertainty. In this study, 6 = 0 is assumed for simplicity. 


III. EXTRACTION OF THE MASS HIERARCHY 


Before the measurement of the surprisingly large sin^ 20 i 3 , it had been known that at the 
oscillation maximum of A 12 , which corresponds to a baseline of approximately 58 km, the 
sensitivity to the mass hierarchy is maximized at sin^ 26*13 ~ 0 . 02 . As sin^ 20 i 3 is no longer 
small, the sensitivity to mass hierarchy needs to be explored as a function of the baseline, 
L, and the detector energy resolution, 6E^/E^. 

In this study, each Monte-Carlo experiment generates a set of 500,000 Ue events by sam¬ 
pling dN°^^/dE°^^ with input parameters, L and dE^/Ei,. The default oscillation parameters 
are taken from Ref. Q fl Q and ..sted .n Table togetben w.h the exp.oted tangea of 
baseline and energy resolution. A total of 72,000 experiment samples are independently 


TABLE I. Default values of neutrino oscillation parameters and the explored ranges of other input 
parameters. 

5 m 2 i[eV^] Am3i[eV^] sin^ 20i2 sin^ 2^13 
7.50 X 10“^ 2.32 X 10“3 0.857 0.1 

L[krri\ a b 

30<L< 60 0<a< 0.08 0 


generated for every 2 km in the baseline and every 0.01 of the energy resolution, dE/^/E^. 
Figure [U shows the Ug L/E spectra at 50-km baseline with the energy resolution varying 
from 0 , which corresponds to an ideal detector, to 0.08/\/E7. As all neutrino masses appear 
in the frequency domain, as indicated by Eq. ( 6 ), a Fourier transform of N{L/Ey) would 
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enhance the sensitivity to the mass hierarchy. The frequency spectrum can be obtained 
using the following Fourier sine transform (FST) and Fourier cosine transform (FCT): 

ftmax 

FST{u) = / N{t) sm{ut)dt, (12) 

^min 

ftmax 

FCTioj) = / N{t) cos{ujt)dt, 

tmin 

where u = 2.54Am?j is the frequency and t = LjE^ is the variable in LjEy space, varying 
from train = L/Emax to tmax = L/Emin- Ouce a huite energy resolution is introduced, the 
phase difference over LjEy is signihcantly smeared out. 


2000 

|S/E = 0 

-SE/E = 0.01/Ie 

-SE/E = 0.02/(e 

1800 

^ jJw 

jpftiL 


1600 

r 

r ta 


1400 


J q 


1200 

- / \ 

r r Vl 

1 H 

1000 


1 ^ 

i A 

800 




600 

1 VUV* ill 

1 Vtfv V 


400 




200 

U . \ 

J . , \ 

J . \ 

2000 

!5E/E = 0.03/(E 

rSE/fE = 0.04/fE 

rSE/E = 0.05/(E 

^ 1800 

- ft. 

a 

A 

^ 1600 


r\ 


(/) 1400 


J \ 

J \ 



- f \ 

f \ 

<u 

J \ 

r A 

J \ 

> 1000 

— f 

“I X 

- r \ 

LU 

J ^ 

f Y — _ 

r \ 

800 

~ I ^ 

~ f X 

f ^ 

O 600 

-I \ 

1 \ 

/ \ 

400 

<0 

1 \ 

/ \ 

/ \ 


— / \ 

- / \ 

- / \ 

cr 

- J , , \ 



2000 

iSE/E = 0.06/fE 

SE/E = 0.07/fE 

SE/E = 0.08/fE 

1800 




1600 




1400 




1200 




1000 




800 




600 




400 




200 


. l/ . . . 1 , . , . 1 . , . , 1 , . , , 1 , , 



5 10 15 20 25 

5 10 15 20 25 

5 10 15 20 25 


L/E [km/MeV] 


FIG. 1. LjEy spectra at 50-km baseline for normal hierarchy (solid line) and inverted hierarchy 
(dotted line) with different detector energy resolutions. 

Figures 2 and 3 show FST and FCT spectra obtained through the Monte-Carlo simulation 
from dm? = 0.002 to 0.028 with energy resolution varied in steps of 2 x 10“^. The impact 
of energy resolution is clear because noisy peaks and valleys fluctuate more with increasing 
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magnitude of energy resolution. The main peak and valley are distinctive and can be used 
to determine the neutrino mass hierarchy while 5Ey/Ey <~ 0.05/\/^. 



0.002 0.0021 0.0022 0.0023 0.0025 0.0026 0.0027 0.0102 0.0021 0.0022 0.0023 0.0025 0.0026 0.0027 0.0i02 0.0021 0.0022 0.0023 0.0025 0.0026 0.0027 0.01 


[eV^] 

FIG. 2. Fourier sine transformed (FST) reactor i>e event rate from 50-km baseline in arbitrary units 
for sin^ 2013 =0.1, normal hierarchy (solid line), and inverted hierarchy (dotted line) at different 
energy resolutions. 


We introduce parameters PVpsT and PVpcT to quantify the features of FST and FCT 
spectra: 


and 


PVpST — 


PVpCT — 


Ap - |A| - 6ml 

Ap + \A^\ \6ml - 6171^1 

Ap — 6ml — 6mp 


(13) 


(14) 


Ap + \ Ay \ \6ml - 6ml\ ’ 
where Ap and A^ are the amplitudes of the peak and valley, respectively, and 5m^ and 5ml 
are the values of 5rn^ at the peak and valley positions, respectively. Figure 4 shows the 
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FIG. 3. Fourier cosine transformed (FCT) reactor Ug event rate from 50-km baseline in arbitrary 
units for sin^ 20 i 3 = 0.1, normal hierarchy (solid line), and inverted hierarchy (dotted line) at 
different energy resolutions. 

distributions of PVpsT and PVpcT for 500 experiments at different energy resolutions. Two 
clusters of points represented by the red open circle (bottom right) and blue solid circle 
(top left) in the plane of {PVpsT, PVpcT — PVfst) corresponding to NH and IH cases, 
respectively, show their own region exactly when SEy/E^ <~ 0.05/\/E^. The upper and 
the lower parts of the scatter plot correspond to IH and NH, respectively. It is shown that 
the distinctive features of NH and IH cases become smeared out as the energy resolution 
worsens, as shown in Fig. 5. 

At large value of sin^ 26^13, the uncertainty of has a little effect on the FST and 

FCT spectra. It comes from the fact that sin^ 2^13 is more effective on narrow modulation 
in the L/E spectrum than sin( 2 |Am 3 ;^|) does in the last term of Eq. ( 6 ). The effect of the 
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uncertainty of is shown in Fig. 6. 
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FIG. 4. PVpcT — PVfst vs. PVfst scatter plots obtained from 50-km baseline for normal 
hierarchy (red open circle) and inverted hierarchy (blue solid circle) at different energy resolutions. 
In the case of an energy resolution of 6E/E < 0.03/^/E, we recognize that points from normal 
hierarchy (bottom right) and points from inverted hierarchy (top left) are well isolated from each 
other. 

Now we consider a method to discriminate between normal and inverted hierarchy using 
the information we gathered from an experiment. We will hnd an experiment on the plane 
of PVpsT versus PVpcT — PVfst as shown in Fig. 4, if we performed the analysis based on 
an approach suggested in this paper. The experiment will be placed on the region of NH 
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FIG. 5. PVpcT — PypsT distributions obtained from 50-km baseline for normal hierarchy (red 
solid line) and inverted hierarchy (blue dotted line) at different energy resolutions. In the case of 
an energy resolution of SEjE > Q.Qh/y/E, we could not distinguish normal hierarchy from inverted 
hierarchy. 


or IH. Could we assess quantitatively whether the mass hierarchy of neutrino has normal or 
inverted hierarchy from the point? This is the probability of being NH (IH) given that an 
experiment happens to be placed on the NH (IH) region: we name it a success probability, 
Pnh{ih)- The probability is simply calculated using classical Bayes’ theorem. For example, 
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FIG. 6. The effect of the uncertainty of on the Fourier spectra at 50-km baseline. Varying 

over its uncertainty range has not significant effect on the characteristic features of Fourier 
sine and cosine spectra. Here lAmg^l = 2.32]']Qgg x 10“^ is considered [l^. 


NH concerned, Bayes’ theorem says, 

r’[x) 

where P{NH\x) is the probability of being NH given that an experiment is found on the NH 
region, P{x\NH) is the probability of being found on the NH reqion given that the hierarchy 
is NH, P{NH) is the probability of being NH, and P{x) is the probability of being found 
on the NH region. 

The probability that an experiment will remain in its own region could be calculated 
from many experiments, which is why we need many experiments. According to classical 
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Bayes’ theorem, there are 


ATNHilH) 

7 ~) success /i 

— nh{IH) ’ 
total 

where is the number of NH (IH) experiments found in the NH (IH) region and 

Ntotai is the number of total experiments found in the NH (IH) region. In this approach, 
50% probability implies a null result. 

Figure 7 shows Pnh and Pjh values obtained from the simulated event samples over 
the baselines at different energy resolutions. We list numerical values of those probabilities 
acquired from MC samples in Table H. In the case of 5E/E < 0.03/\/^, Pnh is greater than 
95% for a baseline of 38-56 km. Similarly, Pjh is greater than 95% at a baseline of 32-52 
km when 5E/E < Q.QA/^/E. As we have no preferred basis to determine which hierarchy is 
correct, we need to introduce a new probability, which shows that an experiment found in a 
region remains in its correct region as long as the energy resolution is sufficient. 

The probability that an experiment will be found in its correct region, namely the deter¬ 
mination probability, Po, is expressed as 


Pd = Pnh ■ Pm- 


( 17 ) 


As shown in Fig. 8 , Pd has a value of > 99% when 5E/E < 0.03/\/^ with a baseline of 
40-52 km. As the energy resolution worsens, Pd rapidly decreases. When 5E/E = 0.04/\/^, 
the baseline is 48-53 km at Pd > 90% as shown in Table H. 


IV. DISCUSSION 

We have studied the experimental requirements to determine neutrino mass hierarchy 
using Fourier sine and cosine transform of the reactor neutrino LjE spectrum at sin^ 26^13 = 
0.1. The parameters PVpsT and PVpcT were dehned to extract features of the Fourier 
sine and cosine spectra, and the mass hierarchy could be obtained from the determination 
probability, Pd based on Bayes’ theorem. 

Since the effect of varying lAmgj^l over its uncertainty has little effect on the FST and 
FCT spectra at sin^ 26*13 = 0.1, the Pd is less dependent on the uncertainty of |Am 3 i| than 
the value of sin^ 2 ^ 13 . 
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TABLE 11. Numerical values of Pnh^ Pih and Pd for O.OS/VE and 0.04/\/E . Errors are calculated 


from binomial distribution. 


baseline [km] 

a 

Pnh 

^Pnh 

Pih 

^PlH 

Pd 

'^Pd 

30 


0.776 

0.024 

0.991 

0.006 

0.769 

0.025 

32 


0.857 

0.020 

0.996 

0.004 

0.853 

0.021 

34 


0.867 

0.020 

1 

0 

0.867 

0.020 

36 


0.938 

0.014 

1 

0 

0.938 

0.014 

38 


0.987 

0.007 

1 

0 

0.987 

0.007 

40 


0.997 

0.003 

1 

0 

0.997 

0.003 

42 

0.03 

1 

0 

1 

0 

1 

0 

44 


0.997 

0.003 

1 

0 

0.997 

0.003 

46 


1 

0 

1 

0 

1 

0 

48 


1 

0 

1 

0 

1 

0 

50 


1 

0 

0.998 

0.002 

0.998 

0.002 

52 


0.997 

0.003 

0.997 

0.003 

0.993 

0.005 

54 


1 

0 

0.935 

0.014 

0.935 

0.014 

55 


0.991 

0.005 

0.907 

0.017 

0.900 

0.018 

56 


0.984 

0.007 

0.846 

0.021 

0.832 

0.022 

30 


0.607 0.0282 0.948 

0.013 

0.576 

0.031 

32 


0.670 

0.027 

0.975 

0.009 

0.653 

0.029 

34 


0.642 

0.028 

0.978 

0.008 

0.629 

0.029 

36 


0.690 

0.027 

0.988 

0.006 

0.682 

0.027 

38 


0.724 

0.026 

0.995 

0.004 

0.720 

0.026 

40 


0.756 

0.025 

1.000 

0.000 

0.756 

0.025 

42 

0.04 

0.867 

0.020 

0.996 

0.004 

0.863 

0.020 

44 


0.847 

0.021 

0.996 

0.004 

0.844 

0.021 

46 


0.861 

0.020 

0.989 

0.006 

0.851 

0.021 

48 


0.905 

0.017 

0.989 

0.006 

0.896 

0.018 

50 


0.913 

0.013 

0.985 

0.005 

0.899 

0.014 

52 


0.938 

0.014 

0.959 

0.011 

0.900 

0.018 

54 


0.891 

0.018 

0.854 

0.020 

0.761 

0.028 

55 


0.900 

P-Pis 

0.811 

0.023 

0.722 

0.029 

56 


0.861 

0.020 

0.816 

0.022 

0.702 

0.030 




FIG. 7. Pnh and Pjh from Eq. (16) as a function of baseline at different energy resolutions. 


As defined in Eq. (16) and Eq. (17), the probability Pd is closely related to the different 
features of each mass hierarchy. Each value of Pd indicates the probability that an experi¬ 
ment will be found inside its correct NH or IH region on the PVfst and PVfct — PVfst 
planes. These different features from different neutrino mass hierarchies suggest that the 
analysis method described in this paper, which is a simple and straightforward approach, 
can be used to determine the neutrino mass hierarchy by using the determination probability 
Pf) based on the Fourier sine and cosine transform of the L/E spectrum. 
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FIG. 8. Determination probabilities, Pd, from Eq. (17) as a function of baseline at different energy 
resolutions. 
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